Bose-Einstein condensation in a stiff TOP trap with adjustable geometry 
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We report on the reahsation of a stiff magnetic trap 
with independently adjustable trap frequencies, lj^ and ujr, 
in the axial and radial directions respectively. This has been 
achieved by applying an axial modulation to a Time-averaged 
Orbiting Potential (TOP) trap. The frequency ratio of the 
trap, iUz/uir, can be decreased continuously from the original 
TOP trap value of 2.83 down to 1.6. We have transferred a 
Bose-Einstein condensate (BEC) into this trap and obtained 
very good agreement between its observed anisotropic expan- 
sion and the hydrodynamic predictions. Our method can be 
extended to obtain a spherical trapping potential, which has 
a geometry of particular theoretical interest. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80. Pj 



Since the first observation of BEC in a dilute alkali 
gas in 1995 0], the spectroscopy of the excitations of 
these novel quantum objects has been an important tool 
for understanding them. A close analogy can be drawn 
between the discrete excited states of a condensate and 
the well defined energy levels of an atom. Spectroscopic 
measurements of atomic energy levels have guided the de- 
velopment of the quantum mechanical model of the atom 
and likewise the spectroscopy of the condensate should 
provide crucial information about its nature. The first 
experiments focused on low energy collective excitations 
both in a pure condensate close to absolute zero Q and 
at finite temperature These latter experiments led 
to a better understanding of the interaction between the 
condensate and the thermal cloud [Q. The first obser- 
vation of a transverse collective excitation was recently 
reported in This also gave evidence for superfluidity 
in a condensate. Recent theoretical work has provided in- 
formation on the behaviour of an excited condensate in 
specific anisotropic trap geometries. This paper presents 
an adjustable magnetic trap that allows us to access these 
geometries whilst retaining a stiff trapping potential. 

In the linear hydrodynamic theory of Bose-Einstein 
condensation, the spectrum of low energy collective ex- 
cited states depends strongly on the trap geometry (i.e. 
the trap stiffnesses in all three orthogonal directions). 
By changing the trap frequency in one specific direc- 
tion, the mode energies change relative to one another 
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and new features appear in the spectrum. For smaller 
condensates, the full Gross-Pitaevskii equation includ- 
ing the kinetic energy term is required. In this regime, 
the size of the condensate, as well as its shape, plays an 
important role in determining the excitation spectrum. 
If the shape of the potential is adjusted so that the sec- 
ond harmonic of one mode corresponds to the fundamen- 
tal frequency of another, then the inherent non-linearity 
of the system leads to coupling between the modes. It 
has been predicted from the hydrodynamic limit of the 
Gross-Pitaevskii equation that this behaviour should be 
observed for a trap potentials with frequency ratios of 
1.95, 1.51 and 0.68 g. 

Another geometry of great interest is the spherically 
symmetric trap, because the excitation spectrum of the 
condensate is simplified by degeneracy and becomes the- 
oretically tractable 0. One mechanism for the decay of 
collective excitations in a condensate is Landau damping, 
which is analogous to a scattering process. A conden- 
sate excitation disappears in a collision with a thermal 
excitation, resulting in a thermal excitation at a higher 
energy. The reduced availability of states in a spherical 
trap, caused by degeneracy, has a major effect on the 
Landau damping rate. Quantitative calculations of this 
process are currently being carried out Q . 

In addition to its importance for studying the collective 
excitations of a BEC, the spherically symmetric harmonic 
trap also produces unique behaviour in a thermal cloud 
of atoms. The breathing mode oscillation is completely 
undamped in this geometry |9j, whilst in an anisotropic 
trap all modes are strongly damped |p^ . 

There are two general types of purely magnetic trap 
for neutral atoms. loffe-Pritchard type traps [0 use 
static magnetic fields to generate a very tight radial con- 
finement. They produce prolate condensates with aspect 
ratios of ~ 0.05. The aspect ratio of a condensate or ther- 
mal cloud is the ratio of the radial to axial size, and is the 
inverse of the trap frequency ratio, in a harmonic poten- 
tial. The TOP trap |12|, in its axially symmetric config- 
uration, combines a quadrupole field with a rotating bias 
field to trap an oblate condensate, with an aspect ratio 
of 2.83. In this paper we demonstrate a novel method 
to independently relax the axial trapping frequency of 
the TOP trap. This new trap, with its adjustable as- 
pect ratio, has been called a zTOP trap because of the 
configuration of bias fields that is used. 

The original attempts to vary the aspect ratio of the 
TOP trap made use of the force of gravity on atoms held 
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in a very weak magnetic potential [Q. This method, 
along with its practical difficulties, is discussed briefly. In 
contrast, the zTOP trap combines an adjustable aspect 
ratio with a stiff trapping potential. The experimental 
setup is described and finally evidence for the successful 
transfer of a condensate into the zTOP trap is presented. 

The standard TOP trap combines a spherical 
quadrupole field of radial gradient B^, with a bias field, 
Bt, rotating in the horizontal xy-plane. The expression 
for the trap potential usually ignores gravity. However, 
for heavy atoms such as ^^Rb, in a weak trap, it must be 
included. The full potential is given by the time-average 
of Eq. 0, where vq = Bt/B'^. 

U{x, y, z, t) = ^lB'q 

( s ^ / N ^ ^1 mgz\ 

\^\ (x + rocoswt<)x + (y + rosmujttjy ~ 2zz\ + j (1) 

At small distances from the trap centre, r ^ rp, this 
approximates to a harmonic trapping potential. The trap 
oscillation frequencies and frequency ratio are 



which has a restoring force in the axial direction that is 
twice as strong as that in the radial direction. Now con- 
sider the effect of applying a TOP bias field rotating at 
ujt in the (radial) xy-plane. This causes the locus of the 
quadrupole field to describe a circle in the xy-plane and 
after time-averaging creates an axially symmetric TOP 
trap. The trap stiffness is reduced in all directions com- 
pared to the quadrupole field trap, but more significantly 
in the radial direction, increasing the aspect ratio to 2.83. 
Extending this argument, the trap stiffness can be pref- 
erentially reduced in the z direction by including an os- 
cillating axial bias field. Our present arrangement uses 
an axial modulation frequency {uja) of 2a;t. Thus in our 
zTOP trap, the locus of the quadrupole field follows a 
three dimensional saddle shape rather than a circle. The 
total zTOP magnetic field is given in Eq. ^ as the sum 
of its three components - the static quadrupole field and 
the oscillating radial and axial bias fields. 
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is the ratio of the gravitational force on the atom to the 
axial magnetic force due to the quadrupole field. Note 
that the trap frequency ratio depends only on rj. Thus 
by measuring the axial and radial trap frequencies, as a 
function of B'^, we were able to accurately calibrate the 
quadrupole field strength. 

However, the aspect ratio only changes significantly in 
very weak traps, where the gravitational force becomes 
important. A ten percent decrease in the aspect ratio 
requires B'^ < 23 G/cm. Since Bt is limited by noise to a 
minimum value of 2 G in our apparatus, this corresponds 
to a radial trapping frequency of < 15 Hz. Attempts were 
made to selectively excite quadrupole oscillations of the 
BEG in this weak gravity trap. However, in such a weak 
trap it was impossible not to excite large dipole motions, 
making the quadrupole oscillation impossible to observe. 

The theory of the zTOP trap is analogous to that for 
the standard TOP trap, except that the bias field os- 
cillates in three rather than two dimensions. Both are 
based on a spherical quadrupole field trap, with a trap- 
ping potential of 



B(<) {xSi + yy-2zz) + 

Bt {cos LUtt X + shlLUtt y) + Bz cos LUat z 



(7) 
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The properties of the zTOP trap are calculated from 
the instantaneous zTOP potential, numerically averaged 
over one cycle of the bias field. This calculation gives 
the following constraints on the frequency, uja, of the ax- 
ial modulation necessary to create a stable trap. Firstly, 
to create a time-averaged potential, LUa must be greater 
than the trap frequencies. Secondly uja must be an ex- 
act integer multiple, n, of the radial (TOP) bias field fre- 
quency, ujt ■ Our calculations showed that deviations from 
this condition produced asymmetry in the time-averaged 
potential between the x and y directions. In the xy- 
plane, the trap will acquire a slightly elliptical, rotating 
cross-section, whilst its centre will describe a small circle. 
These slow variations in the time-averaged potential at 
the beat frequency, - nwj, could heat the cloud. It is 
also possible to choose cot as an integer multiple of uja- 
However this configuration requires the two signals to be 
phase locked, if the 'micromotion' of the trap described 
above is to be avoided. In addition, uja must be lower 
than the Larmor frequency, so that Majorana spin flips 
are not induced. 

The extent of the trap relaxation in the axial direction 
depends on the waveform and amplitude, Bz, of the axial 
modulation, but not on its frequency. The axial modula- 
tion causes the locus of the quadrupole field to oscillate 
in the vertical direction. A waveform that causes it to 
spend more time per cycle at the extremes of this mo- 
tion results in a more relaxed trap in the axial direction. 
Thus an ideal square wave produces the weakest trap, for 
a given maximum Bz field. Removing the higher harmon- 
ics in the axial modulation stiffens the trap as shown in 
Fig.|l|. Note that significant additional relaxation, with 
respect to single frequency modulation, is achieved by 
adding only the third harmonic component of a square 
wave to the fundamental. In this case a spherical trap 
is produced for Bz/Bt ~ 5.5. Since Bt is typically 2G, 
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this is a readily achievable experimental condition. The 
axial trapping frequency is also reduced by increasing the 
modulation amplitude, B^. Thus we are able to vary the 
aspect ratio of the trap during the course of an experi- 
ment by controlling the peak voltage across the axial bias 
coils. 

The apparatus that we use to create ^^Rb condensates 
is described in detail elsewhere In summary, we 

use a differentially pumped double MOT system (which 
incorporates a pyramidal MOT ) to load a TOP trap 
with '--^ 2 X 10* atoms. Evaporative cooling proceeds via 
both Majorana spin flip and radio frequency cutting to 
achieve condensates with ~ 2 x 10'' atoms. 

To generate the axial bias field, two Helmholtz coils 
of 35 turns each, have been added above and below the 
experimental cell. A 2 turn pickup coil is also present 
to monitor the applied field. We currently use an axial 
modulation at 14 kHz (2(jJt). The amplitude of this signal 
is computer controlled and fed to an audio amplifier. A 
transformer is used to match the output of the amplifier 
to the coils. Fine adjustments are made to the frequency 
until the x and z pickup coils produce a clean Lissajous 
figure on the oscilloscope, indicating that the axial oscil- 
lation is at exactly twice the frequency of the radial bias 
field. This frequency relationship is stable to within ^ 1 
Hz over the day. 

The trap frequencies were calibrated with respect to 
the voltage supplied across the axial bias field coils. A 
small, cold thermal cloud was collected in a 'displaced' 
zTOP trap. The additional static magnetic fields causing 
the displacement were suddenly switched off, exciting or- 
thogonal dipole oscillations at the two trap frequencies. 
The radial and axial positions of the cloud were recorded 
as a function of time, using absorption imaging, and fit- 
ted with sine waves to obtain the trap frequencies. Figure 
H shows the data for the normalised trap frequencies and 
the trap frequency ratio, for a range of axial bias field 
amplitudes. 

Theoretical values for the normalised trap frequencies, 
as a function of Bz/ Bt, were found by numerically inte- 
grating the magnetic potential energy of an atom over 
one cycle of the bias field. These values were accurately 
described by an eighth order polynomial. Two free fit- 
ting parameters were used to fit this polynomial to the 
experimental data. One parameter determined the trap 
frequencies without any axial bias field and the other was 
a constant, relating the voltage across the axial bias coils 
to the actual value of Bz/Bt that was produced. Fit- 
ting this expression to the axial trap frequency data gave 
values for both fitting parameters and their respective er- 
rors. The solid lines in Fig.^ show the fitted zTOP trap 
calibration curves. 

The condensate is initially produced by a sequence of 
computer controlled evaporative cooling ramps in a stan- 
dard TOP trap, with final trap frequencies ujr — 126 Hz 
and uJz — 356 Hz (= \/8uJr). The amplitude of the os- 
cillating axial field is then linearly ramped from zero to 
its final value over 0.5s. This adiabatically transfers the 



condensate into the zTOP trap. 

Figure || shows the aspect ratio of a BEG, 12ms after 
release from zTOP traps with a range of frequency ratios. 
Changing the frequency ratio of the trap changes the dis- 
tribution of the self energy of the condensate among the 
three orthogonal directions and hence the form of the 
free expansion of the condensate. The aspect ratio of the 
expanding BEG has been predicted as a function of trap 
geometry using the hydrodynamic equations p6[ | and is 
shown as a solid line in Fig. This theoretical curve con- 
tains no free parameters. The good agreement between 
theory and experiment confirms two facts. Firstly that 
the condensed state survived adiabatic transfer to the 
zTOP trap and secondly that the trap potential experi- 
enced by the condensate has been modified as predicted. 

In conclusion, we have modified our TOP trap to cre- 
ate a stiff zTOP trap with an adjustable geometry. Trap 
frequency ratios between 2.83 and 1.6 have been demon- 
strated by applying a magnetic field in the axial direction, 
oscillating at a single frequency. The addition of higher 
harmonics to the axial bias field modulation will make 
a spherical trap accessible. This prediction, combined 
with the successful transfer of BEGs into the zTOP trap, 
as demonstrated in this paper, promises to open up the 
field of degenerate condensate excitations to experimen- 
tal reasearch. 
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acknowledges the support of a Marie Curie Fellowship, 
TMR program (No. ERB FMBI-CT98-3077). 
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FIG. 1. Simulations of the trap frequency ratio ujz/i^r as a 
function of the Bz/Bt- The solid lines correspond to differ- 
ent Bz waveforms, (A) fundamental only, (B) fundamental -|- 
third harmonic component of a square wave, (C) square wave. 
The dotted line indicates the value of Bz/Bt required for a 
spherically symmetric trap in each case. 
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FIG. 2. Experimental values for the normalised trap fre- 
quencies and aspect ratio as a function of Bz/ Bt, with the the- 
oretical curve that was used to calibrate the zTOP trap, (a) 
shows how the normalised individual trap frequencies, Ur/uiro 
(solid circles) and lj^/l^zO (open circles), depend on Bz/Bt- 
LJrO and LOzo are the radial and axial trap frequencies when 
Bz = 0. (b) shows the aspect ratio of the trap as a function 
of Bz/Bt. For all points, Bt = 2G. 
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FIG. 3. The aspect ratio of a BEG, (Ar/Az), after 12ms 
time-of- flight, as a function of the zTOP trap frequency ratio 
{ujz/cur)- The solid line shows the solution of the hydrody- 
namic equations as a function of trap geometry (no free pa- 
rameters). Each data point is the mean of six shots and the 
error bar shows the standard deviation of the mean. 
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